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FURTHER BUFFETING TESTS IN A CRYOGENIC WIND TUNNEL 

by 


D. G. Mabey (RAE Bedford) 

R. P. Boyden (NASA Langley) 

W. G. Johnson, Jr. (NASA Langley) 


SUMMARY 

Further measurements of buffeting, using wing-root strain gauges, were made in the 
NASA Langley 0.3m Cryogenic Wind Tunnel to refine techniques which will be used in 
larger cryogenic facilities such as the United States National Transonic Facility (NTF) 
and the European Transonic Wind Tunnel (ETW). The questions addressed included the 
relative importance of variations in frequency parameter and Reynolds number, the choice of 
model material (considering both stiffness and damping) and the effects of static aeroelastic 
distortion. 

The main series of tests was made on three half models of slender 65° delta wings with 
a sharp leading edge. The three delta wings had the same planform but widely differing 
bending stiffnesses and frequencies (obtained by varying both the material and the thickness 
of the wings). It was known that the flow on this configuration would be insensitive to 
variations in Reynolds number. 

Additional tests were made on one unswept half- wing of aspect ratio 1.5 with an NPL 
9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic 
speeds. For brevity the test March numbers were restricted to M = 0.21 and 0.35 for the 
delta wings and to M = 0.30 for the unswept wing. 
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1. INTRODUCTION 

The accurate prediction during the design stage of the onset and severity of wing buffeting 
likely to be met in flight is of great importance for both transport and combat aircraft 1 . 
Hence there is considerable interest in developing and refining experimental techniques that 
can be used with confidence during model tests in wind tunnels. 

The dynamic method of measuring wing buffeting on models and aircraft using root-strain 
gauges was suggested for by Huston . Subsequently it was developed by Jones 3 and others 
for use in conventional wind tunnels and in flight. In some preliminary tests, Boyden and 
Johnson 4 showed conclusively that this technique could be employed also in cryogenic wind 
tunnels; also, with such a facility, effects on buffeting due to static aeroelastic distortion and 
variations in Reynolds number could be distinguished. However, some additional tests were 
judged necessary to clarify some outstanding issues. 

Three principal questions were addressed. The first was to confirm that the buffet 
excitation parameter (derived as specified in Ref. 1) could be measured in cryogenic facilities. 
The buffet excitation parameter depends on the total damping (aerodynamic plus structural) 
in the bending mode as well as on the mode shape and generalized mass ( cf discussion of 
equation (3)). Hence this question implies consideration of the variation of both aerodynamic 
and structural damping at cryogenic conditions. Previous experience in a conventional wind 
tunnel suggests 1 that an ideal model would provide a high level of aerodynamic damping, 
and a low level of structural damping because this would facilitate extrapolation of wind- 
tunnel measurements to flight. The present tests imply that appropriate aerodynamic and 
structural levels can be provided on aluminium or steel models in a cryogenic wind tunnel 
(although no steel models were tested). 

The second question relates to the importance of the frequency parameter, n = / c /U . 
For many wings the buffet excitation parameter varies slowly with the frequency parameter 
so that this does not need to be presented precisely. However, for a slender delta wing after 
vortex breakdown, the buffet excitation parameter is sensitive to the frequency parameter 
so that flight full-scale values must be represented precisely. This was confirmed here by 
tests on three half models of slender delta wings with the same planform but with widely 
differing bending stiffnesses and frequencies. The flow on these slender wings was known to 
be insensitive to variations in Reynolds number 6 . 

The third question relates to the suitability of carbon fibre composite materials for 
models to be tested in cryogenic wind tunnels. Models made in carbon fibre have given 
encouraging results for oscillatory tests in conventional wind tunnels because of their stiffness 
and lightness 6,7 . Accordingly two of the delta wings were made in carbon fibre. The 
measurements show that although carbon fibre is a suitable material for models to be tested 
in cryogenic wind tunnels, care is required in laying up the fibres to achieve the required 
bending and torsional stiffnesses. The aerodynamic and structural damping provided by the 
carbon fibre models is satisfactory for buffeting tests. These comparative tests of the delta 
wings also provide information about the effects of static aeroelastic distortion. 

This second series of tests is described in the present paper and definitive recommenda- 
tions are made for future test techniques. This Memorandum stems from a joint research 



project between NASA Langley and RAE Bedford, originally suggested in a lecture at 
Langley 8 . 

2. EXPERIMENTAL DETAILS 

In general the second series of tests had much in common with the original tests 4 . Hence 
for brevity only significant changes are recorded. 

2.1 Wind tunnels and test conditions 

For the original tests 4 the tunnel two-dimensional section was used (breadth = 203 mm 
and height = 610 mm, Fig. la): this section had a slotted floor and roof with solid 
sidewalls. For the present tests a two-dimensional closed adaptive-wall test section was 
used (breadth = 330 mm and height = 330 mm, Fig. lb). The top and bottom adaptive 
walls were kept straight with linear divergence for these tests ( ie the walls were used in a 
nonadaptive mode). This gave a constant Mach number through the test section without a 
model. Only low Mach numbers ( M = 0.35 and 0.21) were used for these tests, which were 
made at total pressures varying between pi = 1.2 bar and 4.8 bar and total temperatures 
varying between 300 K and 110 K. 

2.2 Choice of models 

Two radically different types of half model were tested (Fig. 2). The three 65° delta 
wings (Fig. 2a and Table 1) had a chamfer providing a sharp leading edge. No attempt was 
made to fix transition because of the separations expected from the chamferred leading edge. 
The aluminium delta was the one tested previously 4 . It was 5.1 mm thick and had a first 
bending frequency of 480 Hz. The other two delta wings were made in carbon fibre. One of 
these models had an identical geometry to that of the aluminium delta. However, because 
of the combined effects of the lower density and a small increase in Young’s Modulus, the 
first bending frequency increased to 650 Hz. The other carbon fibre delta has a reduced 
thickness of only 2.5 mm and hence a lower first bending frequency of 360 Hz. Thus with 
these three models, for fixed aerodynamic conditions, both the static bending stiffness and 
the frequency parameter in bending could be varied over wide ranges. 

The bending and torsional frequencies of the three delta models are compared in Table 1 . 
The bending frequencies, taken in conjunction with the wing densities given in Appendix B, 
imply that the bending deflections due to a given static load would be a maximum for the 
thin carbon fibre wing but smaller and about the same for the thick carbon fibre and the 
aluminium wings. This is consistent with the relative bending moduli inferred from the 
deflection measurements, although smaller bending deflections were expected for the thick 
carbon fibre wing (Table 1). In contrast, the torsional frequencies imply that the torsional 
modulus for the thick carbon fibre wing is only 0.74 relative to that of the aluminium wing, 
whereas that for the thin carbon fibre wing is 2.54. The relative static torsional stiffnesses 
inferred were thus 0.74 and 0.30 for the thick and thin carbon fibre wings. However, the 
static twist due to torsion was not measured. Possibly on these delta wings static aeroelastic 
distortion would represent primarily increased wing twist near the tip for both carbon fibre 
composite deltas relative to the aluminium delta. Carbon fibre composite (unlike aluminium 
and steel) is not an isotropic material and differing ratios of torsional and bending stiffness 
ran be achieved depending on how the fibres are aligned. This can be utilized in wing design. 
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The single unswept half-wing (Fig. 2b) was made of aluminium and was tested previously 4 
with free transition. During the present tests transition was fixed by a roughness band of 
micro-beads with a nominal diameter of 0.048 mm extending from x/c = 5% to x/c = 7.5% 
on both the upper and lower surfaces of the wing. The wing first bending frequency was 
270 Hz, as in the previous tests. 

All four models were mounted on the half-model turntable. This was made of the same 
type of the al uminium alloy (7075-T6) as used for the aluminium wings “to minimize the 
effects of thermal expansion and contraction over the range of test temperatures in order 
to provide a model mounting as rigid as possible to keep the structural damping to a 
minimum ” 4 However, use of this turntable precluded the measurement of the overall force 
coefficients, such as lift (Cl), pitching moment (Cm), and drag (Cl>) which could have been 
measured with a half-model balance. These force measurements would have been particularly 
useful when attempting to infer the type of flow separations on the unswept wing. They 
would also have allowed a direct check on the lift on the delta wings given by equation (4). 
In the absence of measurements of the overall lift and normal force on the wings, the wing 
static bending-moment coefficients (derived from the root strain gauges: cf equation (1) in 
section 2.5) become of considerable interest. However, some of these measurements for the 
unswept wing were subject to significant zero errors (cf discussion of Fig. 13). Hence some 
refinement in test technique appears desirable. 

2.3 Instrumentation 

The buffet data acquisition system used for these tests is a two-channel integrated unit 
designed by personnel of the Langley Instrument Research Division. Each channel of the 
buffet system includes switch-selectable gain settings from 200 to 20000, an active low-pass 
filter, a root-mean-square converter, and an integration circuit. It may be operated in either 
an ac or a dc mode. The dc mode is normally used for wind-off calibration of the bending- 
moment strain gauge by application of known moments to the wing with weights. The dc 
mode is also utilized to measure the wind-on, steady-state, root bending moments of the 
wings. In the ac mode, the steady or dc voltages are automatically suppressed. The system 
averages the root-mean-square values of the unsteady voltage signals from the bending- 
moment strain gauge by integrating the signal for a time interval which is preselected from 
1 to 99 s. For this test, an integration time of 20 s was used. The two-pole, low-pass, active 
filter is mounted on plug-in boards and is used to limit the frequency content above the 
range of interest which is the first bending frequency. The unsteady bending-moment signal 
was recorded on a magnetic tape recorder for later off-line analysis. 

2.4 Test procedure 

Before the actual wind-tunnel test, the models were loaded statically while in an 
environmental chamber in order to determine the effect of the temperature range on the 
sensitivity of the root bending-moment strain gauges. The variation in sensitivity with 
temperature was found to be linear for all of the wings over the range from 300 K to 110 K. 
Over this temperature range, the sensitivity for the two aluminium wings increased about 
22% due to a poor match in the type of strain gauges used. For the two carbon-fibre wings 
the sensitivity increased by about 2%. The strain-gauge sensitivity for each model was 
corrected as a function of temperature in the data reduction process. 
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After the wind tunnel had reached the required test conditions and the angle of attack had 
been set, the gain of the buffet measuring system was adjusted to maximize the output of the 
unsteady wing-root bending-moment signal without overloading the instrumentation. This 
signal was monitored with an oscilloscope and a recording oscillograph for amplifier overload 
and for the allowable input range for the analogue tape recorder. The unsteady bending- 
moment signal was then integrated for the 20 s time interval chosen for this test and then the 
integrated voltage and amplifier gain setting were recorded on the tunnel data-acquisition 
system. Afterwards, a 45 s segment of the unsteady signal was recorded on magnetic tape 
for later analysis. Then, the steady root bending-moment signal was measured and recorded 
on the t unn el data-acquisition system. For on-line display purposes, both the steady and 
the dy nami c wing-root bending-moment coefficients were also computed and plotted on a 
desk top micro-computer with an attached plotter. 

The analogue signals recorded on magnetic tape were later digitized at 2500 samples/s 
for off-line analysis. About 43 s of data at each selected test point were processed by digital 
analysis techniques. Measurements of the total damping ratio were made from a least-squares 
fit of the decay of the envelopes of the autocorrelation function^. The frequency of oscillation 
was determined from the average time between peaks of the autocorrelation function. 

2.5 Presentation of results 

The static and dynamic measurements on all four models are presented as coefficients 
defined below. 


Static bending- moment coefficient: 

steady bending moment n ^ 

„o= ’ W 

qb c 

where q is the kinetic pressure, S is the wing area, and c the aerodynamic mean chord. 
Estimates of Cg for the delta and unswept wings are given in Appendix A. 


Cb 


Dynamic bending-moment coefficient: 

rms dynamic bending moment 
qS~c 

The static calibration factor used in equation (1) to obtain steady bending moments from 
steady strain gauge readings is assumed to apply to the conversion of dynamic strain gauge 
readings to the dynamic bending moment. 



Although the coefficient Cq is adequate for a quick assessment of the unsteady measure- 
ments, it is not appropriate to scaling to other structures. This can be achieved only if 
the dynamic measurements are expressed in terms of the buffet excitation parameter in the 
mode excited, defined 1 as: 


\/nG(n) 


2 m 'z 
ypn qS 


c 1/2 . 


( 3 ) 
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where m = generalized mass in mode with respect to motion at tip, 
z — rms of fluctuating tip acceleration in mode, 

£ = total damping — as ratio to critical damping, 
and n — f c /V is the frequency parameter. 

Use of y/nG{n) would have required measurements of the generalized masses and the 
calibration factors between wing-tip accelerations and the wing-root strain signals for all 
four models and these were not available. Estimates of the total damping ratio for the delta 
wings are compared to the measurements in Appendix B. 

3. RESULTS 

3.1 Delta wings 

3.1.1 General description of flow and buffeting characteristics. As a reference, it is 
helpful to recall the main features of the forces and flows (Fig. 3) on 65° delta wings of 
this type observed previously^ in a variable density wind tunnel at ambient temperatures. 
These observations were made on a pair of identical wings (one made of steel and the other 
of magnesium, having widely differing values of Young’s Modulus, E, but the same bending 
frequency) at M = 0.35 and 0.70, in a closed working section selected because of its low level 
of flow unsteadiness. It is important to note that there is no evidence of any effects due to 
static aeroelastic distortion (due to bending or torsion) within these measurements. 

The variation of the normal force coefficient, Cjy, with the angle of incidence, a, is linear 
for both wings from about 0° to 10° (Fig. 3a). The slope, dC N /da, is close to the lift curve 
slope given by slender wing theory of: 


dCL_ r 

da 2 A 


(4) 


From a = 10° to 20° significant non-linear lift develops due to the formation of leading 
edge vortices (eg as shown for a = 20° in Fig. 3b). Fig. 3c shows that the measurements 
of the buffet excitation parameter at the fixed frequency parameter fc/U = 0.34 are the 
gam ft on both wings despite the widely varying values of generalized mass, tip amplitude 
and damping. Formation of the tightly rolled vortices produces a low level of buffeting 
from a = 5° to 15°. Above a = 18° there is a progressive increase in the buffet excitation 
parameter as vortex breakdown moves u pstream and approaches the trailing edge. The light 
and moderate buffeting levels (yj tiG(ti) = 0.0075 and 0.0150) are reached at a — 22 and 26 . 
Ref. 5 does not include any investigation of the effects of frequency parameter on the buffet 
excitation parameter at a constant angle of incidence above vortex breakdown. However 
the measurements of Ref. 4 suggest that above vortex breakdown there is a strong effect of 
the frequency parameter. This is consistent with Keating’s unpublished measurements of 
the excitation spectrum on a slender wing BAC 221 model with vortex breakdown. These 
measurements are important and hence they are reproduced in Fig. 3d from Ref. 1. Although 
the geometry of the BAC 221 was appreciably different from a 65° delta wing, it is known 
that the flows on both wings are broadly comparable. Here the large increase in excitation 
from a frequency parameter, / cq/V , from 0.9 to 1.4 observed in Ref. 4 is consistent with the 
increase in response. We shall see that similar characteristics are found in the present tests. 
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3.1.2 Tests at ambient temperatures. Tests at ambient total temperature at two different 
total pressures on the present three delta wings at M — 0.35 (Fig. 4) provide an interesting 
comparison with the original measurements in a conventional wind tunnel (Ref. 5 and Fig. 3). 
During these more recent tests at constant total temperature, the Reynolds number and 
kinetic pressure increase monotonically with the total pressure, so that for each model the 
effects of any changes due to Reynolds number and aeroelastic distortion are combined, 
as in a conventional variable-density wind tunnel and cannot be separated. (However, the 
subsequent tests will show how the unique capability of the cryogenic tunnel allows the 
separate effects to be identified.) On each particular model the frequency parameter, n, is 
constant because both the frequency and the velocity are constant for this mode of testing. 

For the aluminium wing (Fig. 4a) the static bending-moment coefficients, Cg, are 
identical up to a = 18° for the two pressures, which suggests that the combined effects 
of variations in Reynolds number and static aeroelastic distortion are negligible. When 
account is taken of the displacement of the gauge centre from the model center line (see 
discussion in Appendix A) the slope, dC q / da, is somewhat lower than would be expected 
from the previous tests® (which were consistent with the lift given by equation (4)). For 
angles of incidence above a = 18°, vortex breakdown approaches the trailing edge and the 
static bending-moment coefficients are appreciably lower at the higher total pressure. We 
shall see later (c/discussion of Fig. 5) that this reduction is due primarily to static aeroelastic 
distortion, rather than to the effect of Reynolds number (c/discussion of Fig. 6). The minima 
in Cq (at 20° and 22° at the lower and higher total pressures) are closely associated with 
vortex breakdown. 

For the unsteady measurements the frequency parameter is n = 0.55 and the increases 
in Cb above a = 18° and 20° at the lower and higher total pressures clearly indicate the 
approach of vortex breakdown to the trailing edge. Over the whole incidence range the Cb 
levels are appreciably lower at the higher total pressure. This is because the total damping 
is higher at the higher pressures (£ = 0.025 compared to £ = 0.009) due to the increase in 
aerodynamic damping. Thus the 60% reduction in Cb may be explained by equations (2) 
and (3) under the assumption that there are no static aeroelastic distortion effects and/or 
scale effects on y / nG(n). 

For the thick carbon fibre delta (Fig. 4b) the static bending-moment coefficients are 
identical at both total pressures up to a = 24°. This might be thought to imply that the 
combined effects of static aeroelastic distortion and Reynolds number are negligible over a 
wide incidence range. However the slope dCg/da from a = 0° to 18° is 10% lower than 
for the nominally identical aluminium delta. This large change implies that the lift must be 
reduced, possibly due to the greater twist of the carbon fibre model towards the tip. Possible 
causes are discussed in Appendix C. For this model the static bending- moment coefficient 
above a = 24° is a little higher at the higher total pressure than at the lower total pressure. 
This anomaly is inconsistent with the corresponding measurements on the other two wings 
(cf Figs. 4a and c). For the unsteady measurements the frequency parameter is n = 0.73 
and Cb is again lower at the higher total pressure (because of an increase in total damping 
from C = 0.010 to C = 0.026). 

For the thin carbon-fibre delta (Fig. 4c) the static bending-moment coefficients are 
identical at both total pressures up to a = 8°, with a further reduction in slope compared to 
that of the thick carbon fibre delta. Above a = 8°, the Cb values are appreciably lower at 
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the higher total pressure, consistent with large aeroelastic effects on this flexible model. For 
the unsteady measurements the frequency parameter is n = 0.41 and Cg is again lower at 
the higher total pressure. Note that for this flexible model the C q levels are much smaller 
them for the other two models, again due to the greatly increased aerodynamic damping. 
(The total damping values are £ = 0.020 and 0.060 at the lower and higher kinetic pressures.) 


3.1.3 Tests at constant Reynolds number and Mach number. The unique features of a 
cryogenic wind tunnel were exploited to demonstrate the effects of static aeroelastic distortion 
on the steady measurements, independent of any variations due to Reynolds number. This 
was achieved by tests for a constant Reynolds number of R = 4.8 x 10® at M = 0.35, 
obtained from two different combinations of total temperature and total pressure, giving 
widely different values of kinetic pressure (Fig. 5). On any particular delta, for this mode 
of testing, the frequency parameter is higher at the lower total temperature because of the 
lower velocity. 

For the alumini um delta (Fig. 5a), the static bending-moment coefficients, C#, are 
identical up to a = 18°, which suggests that the effects of static aeroelastic distortion 
are negligible with a small, tightly rolled vortex. The slope, dCg/da, is identical with 
that measured at ambient conditions, which deliberately cover the same range of kinetic 
pressure (c/Fig. 4a). Above a = 18° vortex breakdwon approaches the trailing edge, with 
a clear minim um at a = 22°. Above a = 22° the static bending-moment coefficients are 
appreciably lower at the higher kinetic pressure, indicating a strong influence of aeroelastic 
distortion. This radical change is almost the same as that measured at ambient conditions 
(c/Fig. 4a). This provides an indication that in conventional wind tunnels, effects of static 
aeroelastic distortion could have been confused with, or even attributed to, variations in 
Reynolds n um ber. For the unsteady measurements the signals are small below a = 18°. 
Above a — 18° the buffeting increases as vortex breakdown approaches the trailing edge. 
The much lower values of C g at the higher kinetic pressure are due to the combined effects 
of the change in aerodynamic damping and the change in excitation. The first effect is 
due to the increased aerodynamic damping (proportional to the product pV) which should 
be 2.4 times higher at the higher kinetic pressure. However, this only increases the total 
damping, £, from about £ = 0.020 to 0.025 which would make the response (represented by 
Cjg) about 10% smaller according to equation (3). This is only a relatively small effect. The 
second effect is due to the change in excitation due to the changes in frequency parameter 
and, to a lesser extent, to changes due to static aeroelastic distortion (associated with the 
changed values of C g). For the higher kinetic pressure the frequency parameter is lower and 
this makes the excitation lower. This is consistent with the measurements of the excitation 
spectrum shown in Fig. 3d. 

For the thick carbon fibre delta (Fig. 5b) the static bending-moment coefficients are 
identical up to about a = 12°, which might be thought to imply that the effects of static 
aeroelastic distortion are small over this range. However, the slope, dCg/da, is 10% lower 
than that for the nominally identical aluminium delta. For this model above a — 16° the 
static bending-moment coefficient is appreciably lower at the higher kinetic pressure due 
static aeroelastic distortion, consistent with the measurements on the aluminium delta. The 
unsteady measurements are similar in character to those on the aluminium wing, the Cq 
levels above a = 20° again being much lower at the higher kinetic pressures. Again this 
is due primarily to the lower excitation at the lower frequency parameter. In addition the 
total damping on this model varies from about £ = 0.012 at the lower kinetic pressure to 
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C = 0.025 at the higher kinetic pressure. This increase in damping would make the response 
(represented by Cg) about 30% smaller according to equation (3). This is a significant effect, 
quite apart from that due to the change in frequency parameter. 

For the thin carbon fibre delta (Fig. 5c), the static bending-moment coefficients are 
identical up to a = 8° although with a further reduction in slope compared to that of 
the thick carbon fibre wing. Between a = 8° and 24° the Cb values are appreciably lower at 
the higher kinetic pressure showing a large effect of static aeroelastic distortion, consistent 
with that observed at ambient conditions ( cf Fig. 4c). For a = 28° and 30° the Cb values 
are higher at the higher total pressure. This minor anomaly has not yet been explained and is 
inconsistent with the measurements on the other two wings. For the unsteady measurements, 
Cb is lower again at the higher kinetic pressure, due primarily to the lower excitation at the 
lower frequency parameter. In addition the total damping on this model varies from about 
C = 0.03 to 0.06, which would make the response (represented by Cb) about 30% smaller at 
the higher kinetic pressure. 


3.1.4 Tests at constant kinetic pressure and Mach number. The unique features of a 
cryogenic wind tunnel were exploited to demonstrate the effects of variations in Reynolds 
number at constant kinetic pressure. This can be achieved by testing at a fixed total pressure 
and Mach number (M = 0.35) but varying total temperature (Fig. 6). On a particular 
model, for this mode of testing (provided the elastic stiffness and frequency are invariant 
with temperature), the frequency parameter is higher at the lower total temperature because 
of the lower velocity, as in the investigation of the effects of static aeroelastic distortion (cf 
section 3.1.3 above). The aerodynamic damping is also higher at the lower total temperature, 
as shown previously®. 

For the alu mini um delta (Fig. 6a) there are no scale effects on the static bending-moment 
coefficients over the incidence range from a = 0° to 18°. However, the angle of incidence for 
the minimum in C B increases from a = 20° at R = 1.2 x 10 6 to a = 22° for R = 4.8 x 10 6 . 
This is a precise indication of a genuine scale effect. This effect was observed also in Fig. 4a, 
but there it could not be attributed entirely to the change in Reynolds number because of the 
simultaneous variation in the kinetic pressure, which could have caused aeroelastic distortion. 
For the incidence range a > 24°, there is a further important scale effect. Here the static 
bending-moment coefficients are appreciably higher at the higher Reynolds number. It is 
important to recall that in the corresponding tests at ambient total temperature (Fig. 4a), 
this genuine scale effect was completely hidden by the large opposite effect thought to be 
due to static aeroelastic distortion (Fig. 5a). 

Inevitably the effect of the Reynolds number variation on the unsteady bending-moment 
coefficient is combined with the large effect of the associated change in frequency parameter 
and the small change in damping. Thus according to the unsteady measurements of Fig. 6a, 
the angle of incidence for vortex breakdown increases from about a = 20° at R = 1.2 x 10 6 to 
a = 22° (or even a = 24°) at R = 4.8 x 10 6 . However, the most interesting feature of these 
measurements is the much higher levels at a = 26° and 28° at the high frequency parameter 
associated with the higher Reynolds number. This variation with frequency parameter is 
consistent with the variation in local excitation shown in Fig. 3d. However, some of the much 
smaller increase at the lower Reynolds number is due to the lower level of total damping 
(C = 0.009 at R = 1.2 x 10 6 compared to C = 0.020 at R = 4.8 x 10 6 ). 
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The results for the thick carbon fibre delta (Fig. 6b) show the same general characteristics 
as the al umi nium wing, although differing in significant details. Thus for the steady 
measurements there are no scale effects on the steady bending-moment coefficient over the 
reduced incidence range from a = 0° to 16°. However, it is important to note that the 
slope, dCg/da, is about 10% lower than for the aluminium wing. For this delta there 
is also no minimum in Cg clearly marking vortex breakdown, again probably due to the 
larger distortion. However, above a = 18° the Cg values are appreciably higher at the 
higher Reynolds number, the change suggesting a larger scale effect than that observed 
on the al umi ni um delta. With regard to the unsteady bending-moment coefficient, the 
most interesting feature is the much higher levels of Cg from a = 24° to 28° at the higher 
Reynolds number For this model all of this change can be attributed to the combined effects 
of the variations due to frequency parameter (thought to be large) and the variations due to 
Reynolds n um ber (thought to be small). This is because for this model the total damping did 
not change significantly, being only about C, = 0.011 at both total temperatures. Presumably 
the expected small increase in aerodynamic damping at the lower total temperature is offset 
by a decrease in the structural damping. The absence of any large variation in structural 
damping coefficient with total temperature is an important result (albeit a negative one), 
particularly for the carbon fibre composite models. Further evidence for the absence of large 
variations in structural damping coefficient with total temperature for all models is provided 
in Appendix B. 

All the results for the thin carbon fibre wing (Fig. 6c) are radically different in character 
from those for the other two wings. Thus for the steady bending-moment coefficient there are 
apparently no scale effects over the much wider incidence range from a = 0° to 22°, although 
the slope is appreciably lower than for the thick carbon fibre wing. Above a = 22° there 
is a remarkable change in character, the Cg values now being lower at the higher Reynolds 
number. One tentative explanation for this major anomaly is that at the higher Reynolds 
number the higher lift developed distorts the wing and displaces the vortex inboard. This 
might reduce the static bending-moment coefficient, despite an increase in lift. Whatever the 
explanation, clearly there are major inconsistencies between the steady measurements on this 
model and the other two models. Again with regard to the unsteady measurements there 
are differences in character. The unsteady bending-moment coefficients are low at both 
Reynolds n umb ers, and vortex breakdown is not very clearly marked, even at the higher 
frequency parameter (the higher Reynolds number) where increased excitation would be 
expected. The total damping levels vary from about C = 0.020 at R = 1.2 x 10 6 to £ = 0.032 
at R = 4.8 x 10 6 . Further discussion of these measurements is judged inappropriate. 


3.1.5 Tests at constant kinetic pressure and velocity. With a cryogenic wind tunnel it 
is possible to select a combination of Mach numbers, total pressures and total temperatures 
such that both the kinetic pressure (1/2 pV^) and the velocity (V) remain constant 8 . 
It follows that the product, pV, remains constant, so that for a particular model the 
aerodynamic damping and the frequency parameter will remain constant in the absence of 
Mach number effects. This is an interesting mode of testing for slender wings, where Mach 
number effects are small, at least for attached flows. For slender wings with separated flows, 
the angle of incidence for the formation of the tightly rolled vortices is also independent of 
Mach number, as Fig. 7 shows (based on Fig. 9a of Ref. 4). We have noticed already the large 
changes in the severity of buffeting on a slender wing after vortex breakdown, associated with 
changes in frequency parameter, so that a mode of testing at constant frequency parameter 
has many intrinsic advantages. 
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Fig. 8 illustrates such measurements on the three delta wings at Mach numbers of 
M = 0.21 and 0.35, where the Reynolds numbers are 1.9 xlO 6 and 4.8 x 10 6 respectively. 
For the al umi nium wing (Fig. 8a) the steady bending-moment coefficients are much the same 
from q = 0° to 12°. However between a = 14° and 30° the Cg values at M = 0.35 are 
generally higher than at M = 0.21 due to the combined effects of the increase in Reynolds 
number and Mach number. It should be noted that the differences between the present 
measurements are smaller than those presented in Ref. 4 (shown by the dotted curves in 
Fig. 8a). With regard to the unsteady measurements, vortex breakdown occurs at about 
a = 22° at M = 0.21 and at about a = 24° at M = 0.35. After vortex breakdown the 
unsteady bending-moment coefficients are much the same but somewhat lower than the 
original measurements of Ref. 4 (which were almost identical at the two different Mach 
numbers). For the present tests the total dampings were about £ = 0.012 at M = 0.21 
and about £ = 0.020 at M — 0.35. If the aerodynamic damping were constant (due to 
the constancy of pV ), this increase in total damping must be attributed to an increase in 
structural damping at the lower total temperature. 

For the thick carbon fibre delta (Fig. 8b) the steady bending-moment coefficients have the 
same character but differ significantly, the slope dC g/dct generally being lower at M = 0.21 
than at M = 0.35. In marked contrast, the unsteady bending-moment coefficients are 
virtually identical, with vortex breakdown occurring at about a = 22°. It is significant 
that this excellent correlation of the unsteady measurements occurs when the total damping 
coefficients are almost identical. For this model £ is about 0.016 at M = 0.21 and about 
0.012 at M = 0.35, presumably due to a small decrease in structural damping at the lower 
total temperature if the aerodynamic damping is constant. 

For the thin carbon fibre delta (Fig. 8c), the static bending-moment coefficients are 
identical from a = 0° to 22°, although the slope dCg/da is lower than for the thick 
carbon fibre delta, and much lower than for the aluminium delta. Above a = 22° the 
static bending-moment coefficient is higher at M = 0.21 than at M = 0.35, a result which 
is inconsistent with the measurements on the other models and which appears anomalous. 
However, regarded as a Reynolds number effect, the anomaly would at least be consistent 
with that observed on this wing at M = 0.35 (c/ Fig. 6c). In marked contrast to these 
differences between the steady bending-moment coefficients for a > 22°, the unsteady 
bending-moment coefficients, although small, are virtually identical. Vortex breakdown 
occurs at about a = 24°. Again it is interesting to note that this excellent correlation 
of the uns teady measurements occurs when the total damping values are almost identical. 
For this model £ is about 0.030 at M = 0.21 and about 0.033 at M = 0.35, the small increase 
in total d am ping being attributed to a small increase in structural damping at the lower total 
temperature, on the assumption that the aerodynamic damping remains constant. 

3.2 Unswept wing 

The reader should recall that the present tests were made with fixed transition at 
M = 0.30, to provide a comparison with the previous measurements 4 made with free 
transition over the much wider Mach number range (from M = 0.30 to 0.80). 


3.2.1 General description of flow and buffeting characteristics. No test results are 
available to describe the main features of the separated flow on this wing at low speeds and 
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no flow visualization, or pressure plotting was possible within the present tests. However, 
previous experience on other aerofoils and wings allows a tentative assessment to be made 
(Fig. 9). 

For an NPL 9510 aerofoil at low Reynolds numbers (Fig. 9a), probably a small laminar 
leading edge separation bubble would form at moderate singles of incidence. This would be 
followed almost immediately by turbulent reattachment upstream of the roughness band. 
This tiny bubble would create no significant buffet excitation, or influence the overall forces, 
because of its short length. Thus the first significant excitation would be provided by a 
trailing-edge separation. This would move upstream steadily as the angle of incidence 
increases until the leading edge bubble bursts. With the trailing-edge separation there would 
be a gentle stall and a slow, progressive increase in buffeting. In this low Reynolds number 
range, increasing Reynolds number would reduce the size of the leading edge bubble and 
hence reduce the boundary layer thickness downstream. This would delay flow separation and 
hence reduce (or even eliminate) the buffeting. This would be consistent with the increase in 
the angle of incidence for buffet onset observed on the model wing as the Reynolds number 
increases from R = 0.8 x 10 6 to 3 x 10 6 , with both fixed and free transition (Fig. 10). 

For the aerofoil at high Reynolds numbers (Fig. 9b) the laminar separation bubble would 
burst and extend fairly rapidly downstream, where it would merge with the trailing-edge 
separation moving upstream. If this assessment is correct, the first significant excitation 
would be provided by the leading-edge separation. In this high Reynolds number range the 
mfl.in effect of increasing the Reynolds number is probably to thin the boundary layer in the 
leading-edge region, thus making the effective leading-edge radius somewhat smaller. Hence 
this would tend to provoke leading edge separation and the angle of incidence for buffet onset 
should decrease 1 ®. This would be consistent with the decrease in the angle of incidence for 
buffet onset observed on the model wing as the Reynolds number increases from R = 5 x 10 
to 18 x 10 6 . It is interesting to note that in this region of Reynolds number the angle of 
incidence for buffet onset observed on the model wing is a little higher with fixed transition, 
indicative of a sensitive flow condition in the vicinity of the leading edge. 

Inevitably the three-dimensional flow on a low aspect ratio half wing is more complex than 
on an aerofoil, but it should be similar in character. Fig. 9c sketches the general character of 
the surface flow pattern expected with a trailing-edge separation at low Reynolds number. 
[This sketch is based on observations on a half wing of aspect ratio 1.5 and a NACA 4415 
section given elsewhere 11 .] Fig. 9d shows the general character of the flow expected at 
higher Reynolds number, with a leading-edge separation followed by reattachment and then 
a further separation towards the trailing edge. 


With the separations sketched in Figs. 9c and d, the spanwise variations in the bubble 
length would provide a wide range of excitation frequencies. This excitation when integrated 
in space and time with respect to the first wing bending mode, would be expected to give 
a flat spectrum for y/nG{n) as a function of n. During the previous tests on this wing 
the total damping measurements showed less variation with the product p x V than on the 
delta wing. This was attributed primarily to the higher level of structural damping on the 
unswept wing, together with reduced aerodynamic damping (cf Figs. 13 and 12 in Ref. 4). 
The total damping measurements on the unswept wing in the present tests are not yet 
available but if we assume that the total damping is constant, the values of Cp would be 
expected to be independent of frequency parameter because of the flat excitation spectra. 
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This is as observed in the tests, which will be described shortly. Hence the measurements 
confirm that the excitation on this wing is fairly flat and due to three-dimensional bubble 
type separations, as observed on most other wings of moderate sweep and aspect ratio 1,1 ®. 


3.2.2 Tests at two constant total temperatures. During these tests at two constant tem- 
peratures (Figs. 11 and 12) the Reynolds number and kinetic pressure increase monotonically 
with the total pressure, so that the effects of changes due to Reynolds number and aeroelastic 
distortion are combined. 

For tests at ambient total temperature, T* = 300 K (corresponding with a test in 
a conventional, variable density tunnel) the static bending moment coefficients, C p , axe 
identical at the three Reynolds numbers up to a = 8° (Fig. 11a). This suggests that with 
attached flow on the wing, the combined effects of variations in Reynolds number and static 
aeroelastic distortion are negligible. It is interesting to note that for a < 8° the slope 
dC £ /da is only about 2% less than estimated from aerofoil tests with this section at a 
Reynolds number of about 3 x 10® (See Appendix A and Ref. 11). Above a = 8° there axe 
wide variations with total pressure which are not monotonic with total pressure. We will 
see later that this is because there are significant effects of both static aeroelastic distortion 
(section 3.2.3) and Reynolds number (section 3.2.4). These additional considerations allow 
us to offer an explanation of the curves of Fig. 11a. 

For pt = L2 bar the effects of static aeroelastic distortion are small and there is a 
maximum in C p at a = 16°. For the intermediate pressure [pt = 2.4 bar) the effects of 
static aeroelastic distortion are larger and reduce C p, despite a small favorable scale effect 
due to the increase in Reynolds number. The reduction in Cp is comparable with that 
observed on the delta wing (cf discussion of Fig. 5) and is attributed to a nose down twist of 
the wing outboard sections. For the highest pressure (pt = 4.8 bar) the favorable scale effects 
(due to the reduction of separation) increase Cp greatly and this increase must offset any 
further reduction in Cp due to increased wing twist. These complex changes in Cp suggest 
that in a conventional wind tunnel, operating at ambient total temperature over a range of 
total pressure, it would be virtually impossible to discriminate between the effects due to 
variations in Reynolds number and those due to aeroelastic distortion after flow separation. 

In contrast to the confusion of the steady measurements, the unsteady measurements 
(Figs, llb-d with n = 0.27) show monotonic changes, the incidence for buffet onset being 
respectively 8°, 9° and 13° and pt = 1.2, 2.4 and 4.8 bar. These measurements suggest also 
that for a fixed penetration beyond buffet onset the severity of buffeting is much the same. 
This would be consistent with relatively small changes in total damping, as suggested in 
section 3.2.1. 

For tests at a reduced total temperature, 7) = 110 K, Fig. 12a shows that the static 
bending-moment coefficients are identical at the two total pressures up to a = 11°, which 
suggests once again that on this wing the combined effects of variations in Reynolds number 
and static aeroelastic distortion are negligible. This is consistent with the observation that 
in this range the slope dCp/da is much the same as at T* = 300 K (Fig. 11a) and again 
about 2% less than estimated from aerofoil tests. Above a = 11° there are wide variations 
with total pressure. Here the effects of static aeroelastic distortion (which reduces Cp) and 
the scale effects (which are adverse and also reduce C p) are additive, so that the maximum 
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value of Cjj is much lower at the higher total pressure. The influence of these additive, 
adverse effects on the unsteady measurements (Fig. 12b with n = 0.46) is to lower the 
incidence for buffet onset from a = 13° at pt = 1.2 bar to a = 12° at pt = 4.8 bar. These 
C B measurements also suggest that for a fixed penetration beyond buffet onset the severity 
of buffeting is much the same. This would be consistent with the relatively small changes in 
total damping suggested in section 3.2.1. 


3.2.3 Tests at constant Reynolds number and Mach number. Fig. 13 shows measurements 
for a constant Reynolds number of 3.1 x 10 6 obtained from two different combinations of 
total pressure, giving widely different values of kinetic pressure. With regard to the steady 
bending-moment coefficient, Cg, there is a displacement equivalent to about 2° in the angle 
of incidence over the range from a = 0° to 11° (Fig. 13a). It seems best tentatively to 
attribute this displacement to zero errors in the measurement of Cg because the incidence 
for buffet onset (taken from the_dynamic measurements) is a = 13° at both conditions 
(Fig. 13b). Above a = 13° the C B values are reduced even further at the higher kinetic 
pressure (Fig. 13a), consistent with the trend observed in section 3.2.2. It is remarkable that 
for a fixed penetration beyond buffet onset at a = 13°, the unsteady C B values are much 
the same at both total pressures, despite the changes in frequency parameter, total damping 
and aeroelastic distortion (Fig. 13b). 


3.2.4 Tests at constant kinetic pressure and Mach number. The large influence of 
Reynolds number at low speeds on this wing was illustrated by two tests with constant 
values of kinetic pressure. For the tests at low kinetic pressure (Fig. 14), the total pressure 
was 1.2 bar and the total temperature was reduced from 300 K to 110 K to increase the 
Reynolds number from 0.8 xlO 6 to 3.1 x 10 6 . With regard to the steady bending- moment 
coefficients there are two effects (Fig. 14a). The first is a displacement in the incidence 
range from a = 0° to 8° which is again tentatively attributed to zero errors (as in the 
discussion of Fig. 13a). However above a = 0° much larger values of C B are found at 
the higher Reynolds number, indicative of a strong favorable scale effect which inhibits the 
development of separations. This is consistent with the displacement of the angle of incidence 
for buffet onset from a = 8° to 13°, shown in the unsteady measurements (Fig. 14b). Once 
again, for a fixed penetration beyond buffet onset the severity of the buffeting is much the 

same. 

For the tests at the high kinetic pressure (Fig. 15), the total pressure of 4.8 bar with 
three values of total temperature (300 K, 150 K and 110 K), gives Reynolds numbers of 
3.1 x 10 6 , 7.9 x 10 6 and 12.5 x 10 6 . With regard to the steady bending-moment coefficients 
(Fig. 15a) the displacements of the curves over the incidence range from a = 0° to 12°, are 
consistent with zero errors. However above a = 12° there were large changes which were 
not monotonic with Reynolds number and which must represent serious and complicated 
scale effects. Thus the highest value of C B (at a = 18°) is 0.74 at R — 3.1 x_10 (c/0.84 
at the same Reynolds number but a lower kinetic pressure in Fig. 14a), then C B = 0.76 at 
R = 7.9 x 10 6 but falls to C B = 0.67 at R = 12.5 x 10 6 . Corresponding with these changes in 
the steady bending moment coefficients, for the unsteady measurements there is a systematic 
variation in the angle of incidence for buffet onset. This is a = 13° at R = 3.1 x 10 , increases 
to a = 14° at R = 7.9 x 10 6 and falls to a = 12° at R = 12.5 x 10 6 . It is interesting to 
observe that at a = 18° there are also comparatively large variations in the unsteady values 
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of Cg associated with the variations in the study values of Cg discussed above. Despite the 
difficulties in making unsteady measurements with large separations, these differences in Cg 
are probably genuine and a further indication of scale effect at this high angle of incidence. 


3.2.5 Influence of Mach number. The present tests were restricted to M = 0.30 but 
for completeness brief reference is made to the previous tests, which covered a much wider 
range of Mach number (M = 0.18 to 0.80). The main difficulty with those measurements 
is resolving the disparate effects due to variations in Mach number, Reynolds number and 
static aeroelastic distortion. The most helpful results (Ref. 4, Fig. 10m) are reproduced here 
as Fig. 16, which shows measurements at M = 0.3, 0.6, 0.7 and 0.8 with free transition. For 
M = 0.6, 0.7 and 0.8 the Reynolds numbers and kinetic pressures are not too different and 
hence Mach number effects should predominate over any effects due to aeroelastic distortion 
or Reynolds number. 

The steady bending-moment coefficients (Fig. 16a) show much the same character from 
M = 0.3 to 0.7, the slopes in the attached flow region increasing steadily with Mach number. 
Beyond buffet onset, Cg continues to increase steadily with a. There is a radical difference in 
character at M = 0.80, with a well defined maximum in Cg at 9° (well above buffet onset at 
8°) followed by a minimum at a = 11° (where the buffeting has a maximum). The unsteady 
bending-moment coefficients (Fig. 16b) are much the same in character from M = 0.3 to 
0.7, the angle of incidence for buffet onset falling progressively as observed on many other 
wings 10 . Above buffet onset buffeting increases relatively slowly. However, for M = 0.8 there 
is a radical change in character, which matches that in the steady measurements. Buffet 
onset is sharply defined at a = 8° due to the rapid increase in buffeting. 

Careful examination of the base levels of the Cg measurements shows that they are 
virtually constant over the incidence range from a = 0° to 8°. This is an interesting 
observation because in the incidence range from a = 0° to 4° we might expect trailing-edge 
separation to occur (at least according to the aerofoil tests of Ref. 12, which were made at a 
Reynolds number of about 2.7 x 10 6 ). The absence of such a change in the unsteady response 
must indicate one of three possibilities. The first possibility is that trailing-edge separation 
may not occur at these higher Reynolds numbers. The second possibility is that it does occur, 
but does not alter the wing buffeting. This has been observed by Roos 13 during buffeting 
tests on a supercritical aerofoil at M = 0.82. Such an observation would be contrary to the 
association co mm only believed to hold between trailing-edge pressure divergence (indicative 
of trailing-edge separation) and buffet onset. Here again surface flow visualization would 
have been of interest. The third possibility is that trailing-edge separation does occur, but 
only causes light buffeting which is masked by the tunnel noise. 

4. Discussion 


The measurements presented above are reviewed now in a wider context. Manifestly 
one achievement is the confirmation that in a cryogenic wind tunnel the disparate effects 
of static aeroelastic distortion and Reynolds number can be better isolated whereas in a 
conventional pressurized wind tunnel, operating at ambient total temperatures, these effects 
are inextricably mixed. 
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Thus for the slender delta wings the present tests define the magnitude of scale effects 
and the effects due to aeroelastic distortion at M = 0.35 on both buffet onset and the 
severity of buffeting, together with the large effects of a variation in frequency parameter. 
In addition an indication that the angle of incidence for vortex breakdown increases with 
increasing Mach number can be inferred from the previous tests 4, ^ (cf Fig. 7). Similarly 
for the uns wept wing the present tests define the magnitude of scale effects and effects due 
to aeroelastic distortion at M = 0.30 on both buffet onset and the severity of buffeting, 
together with the small effects of a variation in frequency parameter. An indication that the 
angle of incidence for buffet onset decreases with increasing Mach number can be inferred 
from the previous tests (Fig. 16), consistent with observations on other wings with widely 

varying sections and higher aspect ratios 1 ^. 

When the results are examined closely some additional features are apparent. The 
first feature is that on both types of wing the largest scale effects are in the region of 
buffet onset, consistent with a general conclusion from a review of scale effects in unsteady 
aerodynamics 14 . The second feature relates to the choice of material for dynamic models 
in cryogenic wind tunnels, which lies between aluminium, carbon fibre composite and steel 
(which was not represented in these tests). 

The carbon fibre material withstood the testing at low temperatures and had adequate 
structural damping. However, although the thick carbon fibre model provided a higher 
bending frequency than the aluminium wing, it was not as high as expected. This outcome 
(for the type of excitation spectrum typical of vortex breakdown shown in Fig. 3d), may be 
associated with increased static aeroelastic distortion compared to the aluminium wing. 
Tentatively this distortion has been attributed to twist near the wing tip, which could 
be due to reduced torsional stiffness caused by the lack of isotropy in the carbon fibre 
composite material. This hypothesis is consistent with the static torsional stiffness inferred 
from the wing torsional frequencies (Table 1). It is not clear whether this is peculiar to 
the method of fabricating these particular carbon fibre models. In future applications it is 
possible that carbon fibre models will only be used in cryogenic wind tunnels where high 
frequency parameters are essential {eg when models are being driven, as in the measurements 
of oscillatory pressure described previously^) and when such tests are restricted to low lift 
coefficients. 

By the same reasoning, for general model testing steel alloys would be preferred to 
aluminium alloy because the static aeroelastic distortion will be smaller by a factor of 3, 
deriving from the ratio of Young’s Modulus, E, for steel to that of aluminium. Identical 
models made in steel, aluminium and magnesium have the same frequencies for all modes 
because: 

fa 

and the ratio E/p m is the same for these metals. A steel wing will, of course, give reduced 
aerodynamic damping, considered as a ratio to critical damping compared to an aluminium 
wing but this property (normally regarded as a disadvantage for buffeting tests*) is offset 
partially by the increase in aerodynamic damping for constant pt associated with cryogenic 
temperatures 0 ’ . 

•In contrast, for single degree of freedom flutter the reduced magnitude of the negative aerodynamic damping of the steel wing 
would normally be regarded as an advantage, because this would reduce the amplitude of the limit cycle oscillation encountered. 
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If steel or aluminium is the material selected for models in cryogenic tunnels the frequency 
may be too low on a solid model. This would be particularly serious for slender configurations 
{cf the excitation spectrum of Fig. 3d). However the frequency might be increased by locally 
lightening the model without reducing the stiffness, eg by hollowing out the centre of the 
wing. On the other hand, if the frequency of a steel or an aluminium model were too high, it 
could be lowered by the addition of heavy underwing stores (for a combat aircraft) or by the 
addition of engine nacelles (for a transport aircraft). This principle was discussed in Ref. 8. 

It should be noted that even when made in steel, high aspect ratio model wings will 
suffer significant static aeroelastic distortion near the wing tip, although this will be less 
severe than a full-scale advanced transport aircraft with its flexible structure. Thus static 
aeroelastic distortion will be important at the tip where the local Reynolds numbers are low 
(particularly for winglets) and where accordingly scale effects are likely to be large. The tip 
region often makes an important contribution to the buffet excitation. 

A number of limitations of these tests must be acknowledged. Perhap s the most serious 
limitation is the failure to derive the buffet excitation parameter, y/ nG{n), from equation (3) 
for the al umi nium wings, which had relatively small static aeroelastic distortion. However, 
this would have required measurements of the generalized mass ( m ) and dynamic calibration 
of the strain gauges (giving z). Some measurements of total damping, £, were derived from 
the spectra of the buffeting on the delta wings and these measurements are discussed briefly 
in Appendix B. There is no reason to doubt that if these calibrations had been made, the 
usual levels of buffet excitation parameter observed in wind tunnels 1,16 , and in flight 1,3 , 
would have been achieved. 

Another limitation is the failure to measure the wind-on distortions of the wings and 
their overall forces. Without these measurements the observed anomalies in the C g 
measurements cannot yet be explained {cf Appendix C). (Wind-on distortions have been 
measured successfully in a conventional wind tunnel using a photogrammetric technique 17 
and it should be possible to use this in a cryogenic wind tunnel.) 

Another limitation relates to the increased levels of tunnel noise at the highest frequency 
in the closed two-dimensional working section (Fig. lb) used for the present tests. (This 
potential difficulty was recognized when the tests were suggested 16 .) For all the delta wings 
the low level of buffeting associated with vortex formation (at about a = 5°) is masked 
by the response due to the relatively high level of tunnel noise at that frequency. In the 
original tests 4 on the same aluminium delta wing (at about 490 Hz) in the other (slotted) 
two-dimensional working section, vortex formation also was observed clearly at about a = 5° . 

Fig. 17 compares some typical Cg measurements from the original NASA tests made in 
the top and bottom slotted working section (from Fig. 8d of Ref. 4) with previous RAE 
measurements 5 of /n(?(n ) on the steel and magnesium models taken from Fig. 3c. [The 

scales for yJnG{n) and Cg have been selected to suggest that an affine transformation exists 
between the buffet excitation parameters and the buffeting coefficients. It may be seen 
by inspection of equations (2) and (3) that this would imply virtually constant damping 
for the Cg measurements.) For both series of tests vortex formation is marked clearly. 
For completeness, Fig. 17 includes the levels of yj nG(n) associated with light, moderate 
and heavy levels of buffeting according to well established criteria 1 . Vortex formation is 
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associated with very light levels of buffet at this frequency parameter, and hence cam be 
masked easily by tunnel noise. 

In contrast to the difficulty of determining vortex formation on the delta wings, through- 
out the present tests in the closed working section on the unswept wing (restricted to 
M = 0.3) there was no difficulty in determining buffet onset at / = 270 Hz (the lowest 
bending frequency within these tests). However, in the original tests buffet onset was not 
sharply defined at M = 0.6 and 0.7. This may have implied higher levels of tunnel noise in the 
two-dimensional slotted working section at those speeds, or possibly incipient trailing-edge 
separation. 

Another limitation relates to the uncertainty about the types of flow separation on 
the uns wept wing. When the tests were planned 8 it was recommended that surface flow 
visualization should always be made at buffet onset. However this recommendation could 
not be im plemented for either the present or the original tests^. Although the authors believe 
that Figs. 9c and d give the correct general representation of the type of flow separation on the 
wing, they consider that an attempt should be made to compute the separation boundaries 
on an NPL 9510 aerofoil at M = 0.30. If predictions could be made over the Reynolds 
number range from 1 x 10 6 to say 20 x 10 6 this should suffice to confirm that the separations 
sketched in Figs. 9a and b do occur. 

Despite these limitations, the tests show clearly that the dynamic buffeting technique 
works well in cryogenic facilities. Hence the authors recommend that it should be used 
(as a matter of routine) whenever flow separations are expected on models to be tested in 
cryogenic tunnels. This is particularly important when testing unorthodox configurations 
with unusual planforms (such as an X, W, Z or even diamond) because here the severity 
of buffeting, the magnitude of scale effects and the effects of aeroelastic distortion may be 
radically different from those on orthodox configurations. 

Although half-models offer the advantage of the highest possible Reynolds number, they 
can only give the buffet excitation parameter in symmetric modes. Complete models have 
the disadvantage of providing lower Reynolds numbers (by a factor of l/\/2 for the same 
blockage) but allow the buffet excitation parameter to be determined for both symmetric 
and anti-symmetric modes, as discussed in Ref. 1. In addition complete models are almost 
invariably supported on an internal, six-component strain gauge balance so that overall 
forces should be readily available to complement the buffeting measurements. However, for 
force balance measurements under buffeting conditions, special attention should be paid to 
the model dynamics on the sting so that balance fatigue failures do not occur. Criteria for 
structural damping levels in cryogenic wind tunnels are given in Ref. 19. 

5. CONCLUSIONS AND RECOMMENDATIONS 


This review suggests nine main conclusions and five recommendations with respect to 
buffeting tests in cryogenic wind tunnels. The conclusions are as follows: 

(1) In a cryogenic wind tunnel the disparate effects of variations in Reynolds number and 
static aeroelastic distortion on wing buffeting can be better isolated, whereas in a 
conventional wind tunnel they are inextricably mixed. 
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(2) For delta wings above vortex breakdown the frequency parameter must be correctly 
represented because of the unusual type of excitation. In contrast, for the unswept 
wing the frequency parameter is relatively unimportant. 

(3) Carbon fibre models can be tested successfully in cryogenic wind tunnels and the 
structural damping characteristics are satisfactory. 

(4) Of all the materials (aluminium, carbon fibre and steel) considered for routine wing 
buffeting tests, steel is preferred because of the reduced static aeroelastic distortion at 
the wing tip. 

(5) Carbon fibre may be used when very high frequency parameters must be achieved, eg 
for special tests when models must be oscillated without significant dynamic distortion. 
Careful lay-up of the fibres is required to ensure the requisite bending and torsional 
stiffness. 

(6) For the delta wings tested, with sharp leading edges, scale effects are small up to the 
incidence for vortex breakdown, which itself does not vary strongly with Reynolds 
number for a constant Mach number. 

(7) For the unswept wing scale effects are also small up to buffet onset. 

(8) For the unswept wing buffet onset also varies strongly with Reynolds number in a way 
that suggests that there are two radically different types of incipient flow separation, 
depending upon the range of Reynolds number. 

(9) For all the delta wings and the unswept wing the effects of static aeroelastic distortion 
are small below buffet onset but much larger above buffet onset. However, there are 
anomalous variations in dC^/da for the delta wings which have not yet been explained 
(Appendix C). 

The recommendations are as follows: 

(a) Dynamic buffeting tests should be made as a matter of routine whenever wing models 
are tested in cryogenic wind tunnels. As in a conventional tunnel such buffeting 
measurements help to interpret the steady force measurements, to identify conditions 
suitable for flow visualization, and to identify incipient flow separations on the wing. 
Preferably these buffeting measurements should be expressed in terms of the buffet 
excitation parameter. 

(b) For a given Mach number, these dynamic tests should always include some investiga- 
tion of the effects of variations in Reynolds number, static aeroelastic distortion and 
frequency parameter, to allow precise extrapolation of the measurements to aircraft in 
flight. 

(c) As in conventional tunnels, the levels of flow unsteadiness in cryogenic tunnels should 
be as low as possible and carefully measured. Tunnel operators should ensure that 
conditions ( eg combinations of Mach number and fan speed) which produce high levels 
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of excitation at particular frequencies will not interfere with buffeting measurements 
on the model at the same frequencies. 

(d) When testing unorthodox configurations with unusual planforms (such as an X, W, Z 
or even a diamond) the severity of buffeting, the magnitude of scale effects and the 
effects of aeroelastic distortion may be radically different from those on orthodox 
configurations. Hence these unorthodox configuratio ns shou ld be tested in cryogenic 
tunnels to predict the buffet excitation parameter, y/ nG(n). 

(e) Ideally static bending and torsion stiffnesses should be m easured , perhaps in conjunc- 
tion with the calibrations required to obtain Cg and \/ nG(n). Ideally model static 
and dynamic deflections also should be measured during cryogenic tunnel tests. 
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Appendix A 


ESTIMATION OF STEADY BENDING MOMENTS FOR BOTH WINGS 
AND COMPARISON WITH MEASUREMENTS 

For the slender wings three well-known simple approximations [equations (A-l) to (A-3)] 
were used to estimate the static bending-moment coefficients and thus to identify the effect 
of the outboard displacement of the strain gauges from the root. Corresponding estimates 
for the unswept wing are more difficult, because the bending-moment slope must be inferred 
from the lift-curve slope measured on the same aerofoil. We shall show that these estimates 
are in excellent agreement with the measurements on the unswept wing. 

Delta wings 

According to slender wing theory the load distribution is elliptical and the centre of lift 
is given by: 

y = (A-l) 

However the leading-edge vortices would displace the centre of lift further outboard than 
suggested by equation (A-l). It is reasonable to assume that for the incidence range from 
a = 0° to 10° the normal force coefficient is almost the same as the lift coefficient so that: 

Cn-C L = (^) «• (A-2) 


Hence from equations (1), (A-l) and (A-2): 



(A-3) 


Fig. 4a shows that equation (A-3) gives values much higher than the measurements. 


This discrepancy is influenced by the displacement of the strain gauge centre from the 
model centre line. If the strain gauge centre is taken to be in the middle of the milled 
out area shown in Fig. 2a this displacement is about 0.94 cm or 0.1 s. The effect of this 
displacement is now estimated. 


Suppose the elliptic loading is given by a local lift distribution: 

2Aa 


Cl(v) = — ~y*- 


(A-4) 


Then the normal force outboard of the strain gauge is given by: 

i0=5 


C L = Aa 


8 + 


sin 26 


/ttA\ 

= — <*(0.874), 

sin 1 0.1 V * / 


(A-5) 


rather than by equation (A-2). 
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In addition to the reduction in lift, the moment arm is reduced from (y) to (y — 0.1 s). 
Hence the moment coefficient is given by: 


C B = 


7rA\ 4a s 

7T c 


cos 3 6 


/ sin 20 \ 
_ 0.05 (« + — Jj 


*=5 


0=sin 1 0.1 


( 7 rA\ 4a s 

= (t) sTf 10 - 7791 ’ 


(A-6) 


rather than by equation (A-3). 

Equation (A-6) is closer, but still more than twice as high as the measurements. 


Unswept wing 

A few steady measurements on an NPL 9510 aerofoil in a slotted wind tunnel are available 
at subsonic speeds over a narrow range of angle of incidence where the flow is fully attached. 
From Figs. 4 to 8 of Ref. 12, we find that at M = 0.5 and a Reynolds number of 2.7 xlO 6 
(with transition at 6 to 8% chord on the lower surface and at 4 to 6% chord on the upper 
surface): 

= 6.88/radian. (A-7) 

da 

Then according to the Prandtl-Glauert law for aerofoils: 

J (1 - M 2 ) = Constant. (A-8) 

V da 

The lift curve slope at M = 0 would be lower and given by: 


= 6.88\/l - (0.5) 2 = 5.96/radian. (A-9) 

da v 


If we assume that the interference in the slotted tunnel used for these aerofoil tests was 
small (as the authors of Ref. 11 suggested), equation (A-9) may be used to determine the 
factor /c, which represents the loss of lift on an aerofoil due to the boundary layer thickness, 
according to the relation: 

k = = = 0.95. (A- 10) 


27T da 


6.28 


These results for an aerofoil must now be applied to a wing of finite aspect ratio, A, at 
a Mach number M. For an elliptic wing with an elliptic loading the expression given by 
Glauert for M = 0 can be modified to: 


dC L _ 2vk l (A-ll) 

da Vl - M 2 1 + M/ f_0 


The present half model of aspect ratio 1.5 is equivalent to a complete wing of aspect ratio 3. 
Hence at M — 0.3, from equations (A-10) and (A-ll) we find that: 


^ = 3.76. (A- 12) 

da 
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For an elliptic load distribution, the centre of lift acts as a spanwise distance from the 
centre line given by equation (A-l). 

Hence from equations (A-12) and (1) the estimated value of the bending-moment 
coefficient slope becomes: 

— 3.76 x (0.425) x 1.5 = 2.40/radian. (A-13) 

da da c 

The measured value for the model wing from Fig. 15a for R = 3.1 x 10 6 is about: 

— 2. 37/radian. (A-14) 

da 

Equations (A-13) and (A-14) must be regarded as being in excellent agreement, considering 
the large n um ber of assumptions made in deriving the estimates and the difficulty of making 
the measurements accurately, particularly at cryogenic temperatures. 
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Appendix B 

ESTIMATION OF THE AERODYNAMIC DAMPING FOR THE THREE 
DELTA WINGS AND COMPARISON WITH MEASUREMENTS 


It is convenient to express estimates of aerodynamic damping in terms of the ratio 
to critica) damping. For flat plate wings this introduces a dependence upon the model 
density and thickness. For slender wings of uniform thickness, t, three analytic solutions for 
the aerodynamic damping, 'y, due to vibration in the first bending mode were derived in 
Appendix E of Ref. 5. The preferred expression, from quasi-steady slender wing theory in 
which the aerodynamic damping force on the bending wing is derived from the cyclic changes 


in effective incidence is: 



(B-l) 


where m = model density, 
uj = 27 r/ (radians/s), 

S T = ratio of semi-span/root chord (0.46 for a 65° delta). 

This expression gave good agreement with the aerodynamic damping measured on both 
the steel and magnesium wings in Ref. 5. Hence good agreement would also be expected for 
the present tests. The numerical values needed in equation (B-l) are given in the following 
table. 


Wing 

tu) 

(m/s) 

Density p m 
(g/cm 3 ) 

Aluminum 

15.4 

2.80 

Thick carbon fibre 

20.8 

1.64 

Thin carbon fibre 

5.8 

1.71 


In all three wings the wind-off structural damping is at about the same level, namely: 

g/2 = 0.005. (B-2) 

Since there is no evidence that the structural damping changes significantly due to variations 
in amplitude, normal force or total temperature, the total damping, becomes. 

<= 7+9/2 =(£) (£) °- 903 + 0005 - (B - 3) 
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This theoretical, linear expression is plotted in Fig. 18 together with the measured values. 
The measurements have no clear systematic variation with the angle of incidence (as in the 
tests of Ref. 5) and are scattered about the theoretical expression. The scatter is much higher 
than in the tests of Ref. 5, where the ‘random-decrement’ method was used to determine 
the total damping. In the present tests the measurements may be more scattered due to 
variations in structural damping with amplitude, normal force or total temperature, or to a 
less accurate method of determining the damping, or to a combination of all these effects. 
Variations in structural damping are known to occur — (c/ discussion of Fig. 8a). 

Whatever the cause of the scatter, the aerodynamic damping (which represents the 
dynamic lift due to wing motion in the first bending mode) is adequately predicted by 
slender wing theory. Hence the steady lift should also be accurately described by slender 
wing theory as expressed by equation (4). 
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Appendix C 

SOME UNEXPLAINED ANOMALIES 


Two anomalies have appeared during these tests. The first is that the increase in torsional 
stiffness expected of the carbon fibre deltas was not achieved. The second is that the carbon 
fibre deltas generally give slopes dC q /dct lower than those observed on the aluminum delta. 
These anomalies may be related and the following comments are offered. 

With regard to the loss of stiffness, a value of E in bending about 1.03 x that of aluminum 
was expected from the unidirectional carbon fibre tape used and the lay-up selected (0°, 45°, 
90° and -45°). Table 1 shows that the bending stiffness inferred from the bending frequency 
and the measured densities for the thick carbon fibre was about the same as that for the 
alumi n um wing. This is consistent both with the expected bending stiffness and with the 
values of Young’s Modulus inferred from the static deflections due to load, which were almost 
identical on the al um inum and thick composite delta wings. However the torsional modulus 
inferred from the torsional frequency was less for the thick carbon fibre delta than for the 
aluminum delta (0.76 instead of 1, Table 1). This reduction in torsional rigidity would 
be consistent with larger wing twists at the tip (due to static aeroelastic distortion) which 
would be expected to reduce the C b values. No figure for torsional stiffness is quoted for 
the uni-directional tape. 

With regard to the anomalous slopes, dC B jda , the major problem is in the incidence 
range from a = 0° to 8°. Here the measurements on the three wings generally show’ no 
variation with kinetic pressure, suggesting that the effects of static aeroelastic distortion are 
small. This is a reasonable result, because here the vortex is small and the buffeting low’. 
There are two exceptions when individual wings suggest aeroelastic distortion: the thin wdng 
in Fig. 5c in the range 10° < a < 20° and for the thick wing in Fig. 8b in the range 
0° < a < 18°. 

In contrast to this result for the individual wings, the slope, dCg/da , varies from wing 
to wing in a systematic way, as the table shows (slopes taken from Fig. 5 and the prediction 
according to equation (A-6)). 


Delta Wing 

dC b / do 

Aluminum 

0.40 

Thick CF 

0.32 

Thin CF 

0.30 

Predicted (A-6) 

0.49 


No corresponding changes in normal force curve slope were observed during the tests of 
the steel and magnesium deltas reported in Ref. 5. (Results from those tests are reproduced 
in Fig. 3a). 

It is diffi cult to explain the variation in the present tests. Possibilities considered and 
rejected include: 
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(1) Some subtle model to model difference in the root fixing, which was nominally 
encastre. For the aluminum delta the root block was screwed directly into the aluminum 
turntable. Both carbon fibre wings were bonded to carbon fibre root blocks, tapered in the 
spanwise direction. These were then inserted into a matching taper provided in the aluminum 
turntable. All three end constraints are thought to have provided encastre conditions. 

(2) Some influence of fluctuating bending moments on the measurements of the mean 
bending moment. This must be rejected because the buffeting levels are low in this incidence 
range. In addition a careful test calibration of the instrumentation described in sections 2.3 
and 2.4 shows that fluctuating bending moments do not affect the mean bending moment. 
This assurance is valuable because of the large changes in Cg which occur after vortex 
breakdown (a > 22°). 

(3) Wrong calibration factors for the wings. This explanation must be rejected because 
the static calibration factors have been repeated. The major anomalies occur for the carbon 
fibre deltas and for these the change in calibration factor is only 2% from 300 K to 110 K 
(c/ section 2.4). 

The authors consider that these unexplained anomalies should not be allowed to detract 
from the main achievement: successful buffeting tests in a cryogenic wind tunnel over a very 
wide range of frequency parameter. These anomalies alter neither the conclusions nor the 
recommendations. However these anomalies re-iterate the importance of careful static and 
dynamic calibrations in buffeting tests, and the importance of obtaining good overall force 
data from balances during tests when the effects of variations in Reynolds number and static 
aeroelastic distortion are being evaluated. 


26 



Table 1 


DETAILS OF DELTA WINGS 


Delta wing 

Thickness 

(mm) 

First 

bending 

frequency 

(Hz) 

Relative 
bending 
modulus (£) 
inferred from 
frequency 
and density 

Relative 
bending 
modulus ( E ) 
inferred from 
deflection 
measurements 

Wind-off 
structural 
damping at 
ambient 
conditions 
(fraction 
critical) 

(bending mode) 

First torsion 
frequency 
(Hz) 

Relative 

torsion 

modulus 

inferred 

from 

frequency 

and 

density 

Relative 

static 

torsional 

stiffness 

Thick 

carbon 

fibre 

5.1 

650 ! 

1.07 

1.00 

0.0050 

1630 

0.74 

0.74 

Aluminum 

5.1 

480 

1.00 

1.00 

0.0033 

1450 

1.00 

1.00 

Thin 

carbon 

fibre 

2.5 

360 

1.43 

1.62 

0.0057 

1450 

2.54 

2.54/(2.04) 3 
= 0.30 


Inferred Young’s Modulus for bending for the thick carbon fiber wing is about the same as that for aluminium, te 

E = 73 GN/m 2 (10 x 10 6 lb/in 2 ). 

Specification values for bending of Magnamite Graphite Pre-Preg tape AS4/3502 with unidirectional graphite fibres is 

= 124 GN/m 2 (18 x 10® lb/in 2 ). 

-45° , 0°) which are consistent with the unidirectional values. 




Fibre lay-up orientations are (0°, 45°, 90°, 
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LIST OF SYMBOLS 


A 

Cb 

C b 

Cl, Cb, C m 
Cn 

c 

c 

E 

f 

k 

9 / 2 
M 

771 

^JnG{n) 
n = fc/V 
Pt 

q = 1/2 pV 2 
R 
S 
S t 

s 

T t 

t 

x 

y 


aspect ratio of complete wing 
static bending-moment coefficient (equation (1)) 
dynamic bending-moment coefficient (equation (2)) 
lift, drag and pitching moment coefficients 
normal force coefficient 
local chord 

aerodynamic mean chord of wings 

Young’s Modulus 

frequency (Hz) in bending mode 

factor representing loss of lift on aerofoil due to boundary 
layer (equation (A-10)) 

structural damping coefficient (fraction critical) (equation 

(B-2)) 

Mach number 
generalized mass 

buffet excitation parameter (equation (3)) 
frequency parameter 
total pressure 

free stream kinetic pressure 
Reynolds number based on c 
wing area 

ratio of wing semi-span/root chord (equation (B-l)) 

wing semi-span from centre lin e 

total temperature 

wing thickness 

streamwise co-ordinate 

span wise co-ordinate 
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y 

z 

V 

a 

7 

C 

P 

pm 

U = 27 zf 


distance from centre line of spanwise centre of lift 
rms tip acceleration in mode 
free stream velocity 
wing and incidence 

aerodynamic damping coefficient fraction of critical 
(equation (B-l)) 

total damping, fraction of critical (equation (3)) 
free stream density 
model density 

circular frequency (radians/s) 
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Fig 2 Models used for buffeting tests 
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a) Variation of normal force coefficient with angle 
of incidence 


Fig 3 Main features of forces and flows on delta wings 
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c) Variation of buffet excitation parameter with angle 
of incidence 



d) Spectrum with vortex breakdown 


Fig 3 concluded 











Fig 5 Tests at constant Reynolds number and Mach number (M = 0.35) 







a) Aluminium r/S b) Thick carbon fibre 
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Fig 6 Tests at constant kinetic pressure and Mach number (M 












Tests at constant kinetic pressure and velocity 






Fig 9 
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Fig 9 Tentative assessment of flows near buffet onset at low and high 

Reynolds numbers 















Fig 13 



a) Static b) Dynamic 


Fig 13 Tests at constant Reynolds number and Mach number (M = 0.30) 
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Fig 14 
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a) Static b) Dynamic 


Fig 14 Tests at low kinetic pressure and Mach number (M = 0.30) 
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Fig 17 



Fig 17 Comparison of buffeting measurements on delta wings in two different 

wind tunnels 
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